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Abstract 

In the present paper we introduce some expansions for the Euler 
Gamma function and the Riemann Zeta function. The expansions in- 
volve the falling factorials. In the proofs we use the Faa di Bruno for- 
mula, Bell polynomials, Potential polynomials, Mittag-Leffler polynomi- 
als, derivative polynomials and special numbers (Eulerian numbers and 
Stirling numbers of both kinds). 
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1 Introduction 



Let us first recall some basic facts concerning special numbers and expansions. 
n 

k\ 
partitioning a set of n elements into k nonempty cycles, see [S]). It is set 



By 



we denote the Stirling number of the first kind (number of ways of 



if 71 > 0, 



= 1, 



— for fc > n or fc < 0. The Stirling numbers 



of the first kind fulfill the recurrence formula 

= {n-l) 



n- 1 
k 



n- i 
fc-1 



(1) 



If n > fc + 2 then using formula ([T|) repeatedly, in each step to the last term of 
the resulting sum we get 






k + l- j_ 



(2) 



Stirling numbers of the first kind have the following generating function (see 
Comtet [6], p. 50, p. 135) 



f^ 



(3) 



(i-t)-" = i+ Y. 

l<k<n 

We use common notations for the falling factorial 

{x)k = x{x — 1) ■ ■ ■ {x - k + 1) 
and for the rising factorial (Pochhammer's symbol) 

x^''^ ^x{x + l)---{x + k-l). 

By i?n,fc = B„^k{xi,X2,.-.,Xn~k+i) (see [5] , [5] (p. 133)) we denote the exponen- 
tial partial Bell polynomials in infinite number of variables a;i,a;2,a;3, . . .. We 
recall that the polynomials are defined by the formal double series expansion in 
variables t and u 



exp ( M ^ 



f" 



t" 



^ + Y.—\J2 u''Bn,k{xi,X2, . . .) 



fc=l 



(4) 



By P^ we denote the potential polynomials (see 0, Theorem B p. 141) which 
are defined for each complex number r by 



j-n \ ^ J.71 



(5) 



and 



K ^ Pn(.9i, 92, ■ ■ ■ , Qri) = ^ (r)fcB„,fc(5i,52,...)- (6) 

l<fe<n 

Formulae (O and ^ are a particular case of the Faa di Bruno theorem and P^ is 
the nth derivative (in a point x = a) of the function {G{x)Y ^ where G{x) is given 
as the convergent power series G{x) = 1 + X) 9nV^ /n\, t = x — a, G{a) = 1. 

n>l 

2 The Euler Gamma function 

Substituting in the integral 



r(s+l)= f x^e-'^dx 



X = — log(l — t) we have 



1 1 

r(. + 1) = j{- iog(i - t)ydt = ft^fi+^- + j + 



dt. 



(7) 



Our aim is to find the values of the Bell polynomials Bn.k for the sequence 
(1/2, 2!/3, 3!/4, . . .)■ Using expansion (g]) we get 



exp{uQi+V + V 



= expi ui -1 + 1 + -t+ U^ + U^ + 
= e-" cxpiw (A + \t+ ^t" + -/ + ■■ 



e-^e^^\-{-\og{l-t))\=e-^{l-t)- 



'W E f'S© 



l<k<n '- 



l<fc<n 



1! 2! 3! 4! 



' ^ 1! ^ 2! 3! ^ 4! 



1+ E 
1+ E 



fn-k^k 



7 I 



l<fc<n 



fn-k^k 



E E (-1) 

j=0 l<A;<n 



fn-k^k+j 

n\j\ 



Putting n — k — a>Q,k+j = P>\, {n — k + a, j = (3 — k) we see that the 
coefficient of t^u^ is 



/3 

E 



A: 



(_l)/3-fe 



1 



(fc + a)!(/3-fc)! (a + /?)!^ 



E(-i) 



/3-fc 



fc + a 



a + /3 

;3-fc 



(8) 



Denoting 



Ca,/3 = E(^^- 



p-k 



k=l 



a + k 
k 



a + (3 
13 -k 



we see that the value of the Bell polynomial Ba,fi for the sequence 

(l/2,2!/3,3!/4,...)is 



Ba,P = 



{a + (3)\ 



Ca,l3- 



(9) 



The following table gives the first few values of Ca^p. 

2 3 4 



a\P 1 

1 

10 10 

2 2 3 











3 6 20 15 

4 24 130 210 105 

5 120 924 2380 2520 945 



Lemma 1. Numbers c^^p fulfill the recurrence formula 

Ca,l3 = (a + /3 - l){Ca-l,/3 + Cq_i,^_i). 

Proof Putting into ^ 



P 



Ca-1„ 



a + k — 1 
k 



a + (3-1 
P-k 



k=l 

fc=l 
and adding similar terms we see that formula (jlOp is equivalent to 



3-l-fc 



'a + k -i 
k 



a + (3 -2 
/3-1-fc 



Ca,/3 = (a + /3-l)5](-l)^-'- 

k=l 



a + k-i 
k 



a + (3 -2 
(3-k 



From the other side 



Ca,/3 = X!("l) 



l3-k 



fc=l 

/3 



a + k 
k 



a + 
13 -k 



k=l ^^ ^ 1=1 



a + k — i 
k + 1 -i 



P 



m— 1 



a + (3 — m 
13 + l-m 



E(-ir 



a + l3 



a + /3 - 1 
m — 1 



,_,fa + l3-2 
m — 1 



a + /3- 


- m 


P + l- 


- m 


a^l3- 


- m 


P + l- 


- m 



(10) 



(11) 



(12) 



where we introduced parameter m by m—i = 13— k. Putting a+fc— 1 = a+j3~m, 

(i.e. k ~ 1 — (3 — m) we see that formula p2l) converts to (fTTj) . 

In the start of the above calculation we used formula © and then formulae 

/n-l\ fn-2 

(n — 771 = m — 1 

which are easy to verify. D 

Therefore using ^ and © we see that the integrand in Q can be written 
in the form 



^M 1 + E ^ E(^)^^-^ = ^^ + E ^"^^ E(^)^ 



and finally for the integral ([7]) we get the following expansion 

1 1 s 1 /2 3 , , 

r(s + 1) = \ \ -rS + —s(s - 1) 

^ ^ ' s + 1 s + 2 2 s + 3 V3! 4! ^ ' 

1/6 20 , . 15 , ,, A 

rS + -,s(s-l) + -s(,s-l)(s-2) + 



s + 4 V4! 5! ' '6! 
■J^ + E^EM.T^^ (13) 



Remark 1. Denoting 



.^i" + %tl '^(" + /3)'' 



5a,/3 = (S)^- 



(a + /3)!' 
and using I110\) we see that coefficients ga.p fulfill the recurrence formula 

a + /3-l s- 13 + 1 

9a, = -^— 5q-1,/3 H -^— 5a-l,/3-l- 

a + p a + p 

Thus formula US]) can be rewritten as follows 



r(* + l) = — T + E^E5"./3 (14) 

s + 1 ^^ a + s -f-^ 



a=l ^=1 



3 Derivative polynomials 

Let {ai, 02, . . . , a„} be a permutation of the set {1, 2, . . . , n}. Then {aj, flj+i} is 

an ascent of the permutation if aj < Oj+i. The Eulerian number \ , ) is defined 

as the number of permutations of the set {1,2, ... ,n} having k permutation 
ascents (see [5], p. 267). For example for n = 3 the permutation {1, 2, 3} has two 
ascents, namely {1,2} and {2, 3}, and {3, 2,1} has no ascents. Each of the other 

four permutations of the set has exactly one ascent. Thus (|-,) = 1:( )=4, 

Consider a function x = x{t) which satisfies Riccati's differential equation with 
constant coefficients 

x' (t) ^ ax'^ + bx + c, (15) 

where a, b, c are real numbers, a ^ 0. Examples of such functions and equations 
are: 

1. x{t) = tani, x'{t) = x'^ + 1, 

2. x{t) = tanhi, x'{t) = -x^ + 1, 

3. x(i) = l/(l + e*), x'{t)^x^-x. 

Let ax^ + bx + c — a{x — a){x — /3), where the roots a, j3 are real or complex 
conjugate numbers. In paper |14j the following theorem is proved. 

Theorem 1. // a function x{t) satisfies equation i tJ5|J 

x' {t) = a{x - a){x - P) (16) 

then the nth derivative of x{t) can be expressed by the following formula: 

n-l , > 

x(")(t) = a" ^ ^^^(x - a)^-+i(x - /?)"-^ (17) 

where n — 2,3, . . .. 

The polynomials, of the variable x, on the right hand side of (|17p are a partic- 
ular case of the so called derivative polynomials which were recently intensively 
studied (see for example [1], [5], [7], [5], [Hj, [13] )• The derivative polynomial 
of the nth derivative of the function x{t) = 1/(1 + e*) (in ([T5|) a = 0, /3 = 1), of 
degree n + 1 will be denoted by Qn+i{x). Thus for n > 2 

n-2 , _, \ 

'n — 1\ 



fc=0 ^ ' 



It is seen that, for any n > 2, the polynomial Qn{x) is divisible by x{x — 1). 
The polynomial, of degree n — 2, resulting from the operation will be denoted 
by Pn-2(x). Thus for any n > we have 



p„(.) = ^(^" + ^y(x-ir 



-k 



In paper fTSJ (see also [1], [S], [7]) it is proved that the polynomials Q„(x) can 

be expressed in terms of the Stirling numbers of the second kind I > (number 

of ways of partitioning a set of n elements into k nonempty subsets, see Graham 
et al. [5]), namely 



k=l ^ ^ 



Then one can easily check that 

n+l 



p„(.) = 5:(-iri-'=fc!|"+^| 



x^-\ 



By results of [13] it follows that for n = 1, 2, . . . the polynomial P„ [x) has exactly 
n simple zeroes which all are real and lie in the interval (0, 1). Moreover the 
zeroes of the polynomials have the interlacing property. 

4 The Riemann Zeta function 

It is well known that if 5fts > then 



-1^ 

1 + e* 



t'-^^—^dt = C{s){l-'2'-')V{s). (18) 





Integrating n + l times by parts the left hand side of p^ we obtain successively 

OO CO oo 

t'-'-^dt ^ ~- ft^ (t^) dt ^ ^^ /i^+1 f — L-) dt 

1 + e* sj Vl + e*/ s(s + l)7 Vl + eV 



= ■■■= , \,^ , 7/*'+" T^ '^^- (19) 

s(s + l)---(s + n) 7 Vl + eV 



Substituting in (fT9)) . x — 1/(1 + cxp(i)) we have 

s(s + 1) • • • (s + n) 7 \ X / x[x-\) 





1/2 

= f ^u ( ^ \ / log ^"(^) '^^' (20) 

s(s + 1) ■ ■ • (s + n) J \ X J 



where the polynomials Qn{x) and Pn{x) have been defined in the previous 
section. 

Our next aim is to find a power series expansion of the function (log((l — a;)/a;))'' 
by using the Potential polynomials ^. Since 

>°^ ^ - '»« 1^1^ - «> - ^^•> (' + 5" - ^"^'' 4» - ^■^'' + ■ ■ > 

we have to find the expansion of the function 

1,1 

In order to do it we will compute the Bell exponential polynomials 
B„,fc(0, 2!/3, 0, 4!/5, 0, 6!/7, . . .)• By using © we get 

1- 1.6 



expL Qt^ + -f" + -e + 

= e"" expj u M + it^ + it* + it^ + 

= e-"exp{^2(t+V + V + V 

f u , l + t1 _,, /1 + t 

= e exp^; —log 



^2t 1-tJ VI -^ 

It is well known that the function ((1 + t)/{l — t))^ is a generating function for 
the Mittag-Leffler polynomials Mfc(x) 

^ fc=0 

A first few Mittag-Leffler polynomials are as follows 

Mo{x) = 1, 

Mi{x) = 2x, 

M2{x) = Ax"", 

Msix) = 8a;^ + 4a;, 

M4(x) = 16a;'' + 32x'^, 

M^{x) = Z2x^ + 160a;3 + 48a;. 

Bateman [1] regards polynomials gk{x) = Mk{x)/k\ and gives for them the 
following recurrence formula 

ngn[x) = {n- 2)gn-2{x) + 2xgn-i{x), (21) 



with which he refers to Bclorizky [3]. Multiplying both sides of (PT|) by (n — 1)! 
we get 

A/„(x) = (n - l)(n - 2)M„_2(a;) + 2a;Af„_i(a;). (22) 

n 

If Mn{x) = J2 o.n,kx'' then the recurrence formula ([^ yields 

fe=0 

fln.fc = (n - l)(n - 2)a„_2,fe + 2a„_i,fe_i. (23) 

The following table gives the first few values of a„_fc. 

6 7 



n\k 





1 


2 


3 


4 


5 


6 





1 















1 





2 













2 








4 











3 





4 





8 









4 








32 





16 







5 





48 





160 





32 





6 








736 





640 





64 



7 1440 6272 2240 128 

Lemma 2. Numbers a„^fe fulfill the following formula 

fe-i 
a„,fc = ^ 2™(n - 1 - m)(n - 2 - m)a„_2-m,fc-m- (24) 

Proof. Formula ([M|) follows by using formula (|23p repeatedly, in each step to 
the last term of the resulting sum. D 



Therefore 

OO -, 71 



l—t) -' ^ n\ ' "' Z.;^2.«».H2^J 

J=0 -^ ■ y n=0 fe=0 






"" ' ,n-k 
2,k 2fc '' 



oo a P -. -. ^ 

Q=0/3=0 ^ ^^ j"=0 \ J / 

and we see that, in our case, the value of the Bell polynomial Ba^p is 



Let us denote 
















/3 
J=0 


\ J ) 


1 




The following table 


gives the first few values of b^^p. 






a\(3 




1 


2 


3 


4 


5 







1 










1 















2 




2 











3 
















4 




24 


40 









5 

















6 




720 


2688 


2240 







7 


















8 




40320 


245376 


443520 


246400 






(26) 



Lemma 3. Numbers ha^p fulfill the recurrence formula 

bc.,p = {a + l3~ 2){a + l3 - l)(6„-2,^ + &a-2,/?-i). (27) 



Proof. Assume that a > 4. Since the last term in (p6)) is zero then as the upper 
limit of the sum we can take (3 — 1. We have 

X-^\ ,,/a + ^~2\ 1 , , 

f'a-2,/3 = 2^(-l)-'l j^p— a»+fl-2-j,0-j, (28) 

13-2 

By adding similar terms in (f28| and ([29]) we see that formula (f27)) is equivalent 
to 

&a,/^-(a + /3-2)(a + /?-l)^(-l)^"(" ''' )__a„+^_2_^.^_^, (30) 



7 v^ / 1 M- / ct + /5 ~ 3\ 1 ,„„, 

&a-2,/3-l = 2^ (-!)■'( I 0/3-1-J- «c.+^~3-3-,ff-3-- (29) 



^ J 



From the other side 

J = \ ^ / j^Q \ J / 

X ^ 2''{a + 13 - j -i-l){a + li - j -i-2)aa+p-j-.i^2,p-3- 

i=0 



10 



/3-1 



X! oF^*^" + ^ - m - l)(a + ^ - m - 2)aa+0_,„_2,/3-m 



m=0 

T7l 



'^^^ f— 1)"' /a+B—l\ 



X 
k=0 



m=0 



= (a+/3-l)(a+/3-2)^ 



(-1)™ fa + l3-3 



m=0 



aa+/3-m-2,/3-m, 



where tti = i + j- Thus formula ([27]) is proved. In the start of the above 
calculation we used formula (p4|) and then the following formulae 



)(n-m- l)(n-m-2) = ("^ |(n-l)(n-2), 
TO / \ m J 

which are easy to check. D 

Using formula ([5]) and the values of the Bell polynomials (|25p we get 

l-xV' / l + (l-2a; 
log = log 



2 fl - 2x) 



1 - (1 - 2x) 
(l-2a;)3 (l~2xf 

^2'xi-2.r(i+(i^ + (l^ 

= 2'-(l-2xr|l + |(l-22;fr+(l-2a;)4(^^r + ^r(r-l) 

/, „ xB /^720 2688 , , 2240 , ,,, „A 
+ (1 - 2xf (^— r + ^^r{r - 1) + -^r{r - l)(r - 2) j + ■ ■ • 

/ oo L™/2J , 

.r(i-2.r i+i:(i-2.,-i:^(.-). . 

^ n=l /£=1 ^ ' ' 

Therefore since 

1/2 

(1 - 2x^1 dx 



2(7+1) 





11 



then putting in (I^Ul) say n ~ {Pq{x) = 1) we get 
1 /720 2688 , , 2240 , 

,s_l / -, <x) _, L™/2J 



^ s + m+1 ^ (m + fc)!^''^^'' 

m=l fe=l 

TT ^ s + 2to + 1 f^ (2m + A:)! '^'^ 



s \s + l ^ s + 2m+l^^ (2m + k)\ 

^ m=l k=l ^ ' 

Remark 2. Denoting 

Zm,k = 7 ; T\\\^l^-' 

[m + ky. 

and using the recurrence formula for bm,k we see that coefficients z^^k fulfill the 
recurrence formula 

m + k — 2 s — fc + 1 

Zm,k = —, Zm~2,k H —, Zm-2,k~l- 

m + k m + k 

Thus the above expansion can be rewritten as follows 

2S-1 



C(.)(i-2i-)r(s) = 



m—1 k—1 



Similar expansion can be established for any n > 1 in (1201) . In fact the 
expansion can be seen as an expansion for the integral 



(^log— ^j P„(x)dx = ^^^C(s)(l-2i-)r(s). (31) 



If we fix an s £ C such that 1/2 < 5Rs < 1 and ((s) ^ then it is seen 
from (f3T|) that the left hand side (i.e. our expansion of it) goes to infinity for 
n — >■ oo. Therefore if for any s G C such that 1/2 < 5fts < 1 we had indicated 
a nonnegative integer n such the evaluation, at that s, of the expansion of the 
left hand side of (|3T]) is different from zero then the RH would follow. 
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